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Abstract. We illustrate the use of the statistical method of moments for determining the 
position and momentum distributions of a quantum object from the statistics of a single mea- 
surement. The method is used for three different, though related, models; the sequential 
measurement model, the Arthurs-Kelly model and the eight-port homodyne detection model. 
In each case, the method of moments gives the position and momentum distribution for a large 
class of initial states, the relevant condition being the exponential boundedness of the distribu- 
tions. 

PACS numbers: 03.65.-w, 03.67.-a 



1. Introduction 

One of the main problems of quantum mechanics deals with the possibility of measuring 
together the position and momentum distributions and of a quantum system prepared in a 
state p. The basic structures of quantum mechanics dictate that there is no (joint) measurement 
which would directly give both the position and momentum distributions and that, for instance, 
any determination of the position distribution necessarily disturbs the system such that the 
initial momentum distribution p^ gets drastically changed. 

In recent years two important steps have been taken in solving this problem. First of all, 
the original ideas of Heisenberg [lO] have finally been brought to a successful end with the 
seminal paper of Werner |33j which gives operationally feasible necessary and sufficient con- 
ditions for a measurement to serve as an approximate joint measurement of the position and 
momentum distributions, including also the inaccuracy-disturbance aspect of the problem. The 
second breakthrough in studying this question comes from a reconstruction of the state p from 
a single informationally complete measurement, notably realized optically by an eight-port 
homydyne detection [20], |2T[ p. 147-155] (for a rigorous quantum mechanical treatment, see 
[T5]). In conjunction with an explicit state reconstruction formula (known at least for the 
Husimi-distribution [7]), this allows one to immediately determine the distributions of any 
given observables. 

If one is only interested in determining the position and momentum distributions and p"^, 
it is obviously unnecessary to reconstruct the entire state; one should be able to do this with 
less information. Here we will use the statistical method of moments to achieve a scheme for 
position and momentum tomography, i.e. the reconstruction of the position and momentum 
distributions from the measured statistics. The price for using moments is, of course, that they 
do not exist for all states, and even when they do, they typically do not determine the distri- 
bution uniquely. Hence, we restrict here to the states for which the position and momentum 
distributions are exponentially bounded. We note that this is an operational condition and 
can, in principle, be tested for a given moment sequence |19| . 

We consider three different, though related, measurement schemes based on the von Neumann 
model, Sect. 13.11 and the balanced homodyne detection technic, Sect. 14.31 The first model 
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is a sequential measurement of a standard position measurement of the von Neumann type 
[30] followed by any momentum measurement, Sect. 14.11 The second (Sect. 14. 2p builds on 
the Arthurs-Kelly model [2] as developed further by Busch [4] whereas the third (Sect. 14. 3p 
model uses the quantum optical realizations of position and momentum as the corresponding 
quadrature observables of a (single mode) signal field implemented by balanced homodyne 
detection |T4]. In Sect. [5] we apply the method of moments to determine both the position and 
momentum distributions and from the actually measured statistics. Finally, we compare 
our method with the state reconstuction method. Sect. [H There we also comment briefly the 
possibility of inverting convolutions. We begin, however, with quoting the basic no-go results 
on the position-momentum joint /sequential measurements. 

2. No JOINT MEASUREMENTS 

There are many formulations of the basic fact that position and momentum of a quantum 
object cannot be measured jointly, or, equivalently, that, say, any position measurement 'de- 
stroys' all the information on the momentum prior to the measurement. In this section we 
recall one of the most striking formulations of this fact. To do that we fix first some notations. 

Let 7i be a complex separable Hilbert space and L{7i) the set of bounded operators on H.. 
Let be a nonempty set and A a cr-algebra of subsets of VL. The set function E : A ^ L{7i) 
is a semispectral measure, or normalized positive operator measure, POM, for short, if the set 
function A 3 X {ip\E{X)'il}) G C is a probability measure for each ip & Hi, the set of unit 
vectors oiTC. We denote this probability measure by p^. A semispectral measure E is a spectral 
measure if it is projection valued, that is, all the operators E{X),X G A, are projections. IfTC 
is the Hilbert space of a quantum system, then the observables of the system are represented by 
semispectral measures E and the numbers {ilj\E{X)ip) , X ^ A, ip & Hi, are the measurement 
outcome probabilities for in a vector state ip- An observable is called sharp if it is represented 
by a spectral measure. Otherwise, we call it unsharp. Here we consider only the cases where 
the measurement outcomes are real numbers, that is, {Q,A) is the real Borel space (M, i3(]R)), 
or, pairs of real numbers, in which case {fl,A) is (M^,;B(M^)). The position and momentum 
distributions p'^ and p^ are just the probability measures and p^ defined by Q and P together 
with a density matrix (mixed state) p. 

An observable M : i3(M^) L(H) has two marginal observables Mi and M2 defined by the 
conditions Ml (X) = M(XxM) andM2(y) = M{RxY) for allX,F G B{R). Any measurement 
of M constitutes a joint measurement of Mi and M2. On the other hand, any two observables 
El and E2 admit a joint measurement (or equivalently a sequential joint measurement) if there 
is an observable (on the product value space) M : B(M?) L(H) such that Ei = Mi and 
E2 = M2. The following result is crucialu 

Lemma 1. Let M : B{M.'^) L(H) he a semispectral measure, such that one of the marginals 
IS a spectral measure. Then, for any X,Y e BiM), Mi{X)M2{Y) = M2{Y)Mi{X), that is, 
the marginals commute with each other, and M{X x Y) = Mi{X)M2{Y) , that is, M is of the 
product form. 

Assume that M : B(M^) — > L{H) is an observable with, say, the first marginal observable 
Ml being the position of the object. Then Mi and M2 commute with each other, and due 

^This result seems to be well-known, and part of the proof goes back to Ludwig [221 Theorem 1.3.1]. However, 
we were unable to identify a full proof in the literature, and so we give one in the appendix 
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to the maximality of the position observable Q any M2(Y) is a function of Q. Therefore, M2 
cannot represent (any nontrivial version of) the momentum observable. Similarly, if one of 
the marginal observables is the momentum observable, then the two marginal observables are 
pairwisely commutative, and the effects of the other marginal observable are functions of the 
momentum observable. 

3. Position/momentum measurements 

It is a basic result of the quantum theory of measurement that each observable (sharp or 
unharp) admits a realization in terms of a measurement scheme, that is, each observable has 
a measurement dilation [26]. In particular, this is true for the position and momentum ob- 
servables Q and P. However, due to the continuity of these observables they do not admit 
any repeatable measurements [26l [23]. In fact, the known realistic models for position and 
momentum measurements serve only as their approximative measurements which constitute 
Q and P -measurements only in some appropriate limits. Here we consider two such models, 
the standard von Neumann model and the optical version of a Q, resp. P, -measurement in 
terms of a balanced homodyne detection. Before entering these models we briefly recall the 
notion of intrinsic noise of an observable and the corresponding characterization of noiseless 
measurements. 

For an observable E : i3(]R) — L{H) the k^^ moment operator is the (weakly defined) 
symmetric operator E[k] = J-^x^ dE with its natural (maximal) domain D[E[k]). In particular, 
the number {il)\E[k]'il)) = J^x^ dp^ is the fc**^ moment of the probability measure p^. The 
(intrinsic) noise of E is defined as N{E) = E[2] — and it is known to be positive, that 

is, {ip\N{E)<^) > for all <^ e D{E[2]) n D{E[lf). If the first moment operator E[l] of E is 
selfadjoint, then E is sharp exactly when E is noiseless, that is, N{E) = [16j@ We recall also 
that the first moment operator E[l] of an observable alone is never sufficient to determine the 
actual observable. In statistical terms, the first moment information (expectation) {ip\E[l]ilj) , 
ip eHi, does not suffice to determine the measured observable E. 

3.1. The von Neumann model. Consider the von Neumann model of a position measure- 
ment of an object confined to move in one spatial dimension [301 Sect. VI. 3], see also e.g. [5], 
Sect. II. 3. 4]. Let TC = L^(M) be the Hilbert space of the object system, and let Q denote its 
position operator. We let Q denote the spectral measure of Q. To measure Q we couple it with 
the momentum Pq of the probe system, with the Hilbert space /C = i^^(M), and we monitor the 
shifts in probe's position Qq, with the spectral measure Qq. Let U = e~*'*''3®^o unitary 
measurement coupling, with a coupling constant A>O,0G/C, ||0||=1, the initial probe 
state, and let : ?i — > ?i ® /C denote the embedding ^^(v?) = if ® (j). The actually measured 
observable of the object system is then given by measurement dilation formula 

E{X) = V*U*I ® Qo(X)f/r^, X e S(R). 



The selfadjointness of the first moment operator is crucial for this condition. Indeed, if, for instance, one 
restricts the spectral measure of the momentum observable P in L^(R) by a projection Q(/),/ = [a,b], to get 
a POM P -.Y ^ Q(/)P(r)Q(/)|i2(/) acting on L'^{I), one has P[fc] = P[l]'^ for all k, and thus also 7V(P) = 0, 
though the first moment P[l] is only a densely defined symmetric operator [8]. This is also an example of the 
variance free observables as discussed in [32] . A noiseless observable is variance free, but due to the domain 
conditions the reverse implication may not be true. 
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A direct computation shows that E is an unsharp position, with the effects 
(1) E{X) = {xx*f){Q), 

where Xx * f denotes the convolution of the characteristic function xx of the set X G B{M.) 
with the probability density /(x) = A|0(— Aa;)p. 

3.1.1. Limiting observable. The actually measured observable E depends on two parameters: 
the coupling constant A and the initial probe state 0, that is, E = E^''^'. The structure of 
the effects jl]) suggests that the semispectral measure E comes close to the spectral measure 
Q whenever the convolution xx * f comes close to xx- This evident fact can be quantified in 
various ways. 

Due to the convolution structure of E, the geometric distance between the observables E 
and Q can easily be computed [33], and one finds that 

d{E,Q) = jj \x\\(l){x)\^dx, 

showing that whenever the integral is finite, the geometric distance tends to zero as A increases, 
or becomes more sharply concentrated around the origin. It follows from the definition 

of the geometric distance, that d{E, Q) = implies E = Q. However, this does not settle the 
question of the limit -E Q in either of the two possible intuitive meanings. For that we use 
the method of moments. 

In order to be able to determine the moment operators of the unsharp position observable E, 
we assume that G Cj^(M), so that, in particular G D{Qq) for each G N. In that case the 
moment operators E[k] can all be computed]^ and they turn out to be polynomials of degree k 
of g, that is, D{E[k]) = D{Q^), and 



(2) m = t (*) 



Therefore, in particular, on D{E[2]) = D{Q^), one has N{E) = p-Var (Qo, 0)-^, suggesting, 
again, that, for a fixed 0, if A is large, then the noise N{E) is small, or, for a fixed A, if 
Var {Qq, 0) is small, then, again, N{E) would be small. But, again, the precise meaning of the 
limit E ^ Q in either of the cases A — > oo or Var {Qq, 0) — > waits to be qualified. 

Consider first the limit A — oo, so that, the operator measures are actually E^, with the 
moment operators E^[k] of (l2|). Let D be the linear hull of the Hermite functions, so that 
D C D{Q'') = D{E^[k]) for all k (and for all A), and 

(3) lim {ij\E'[k]^) = (V^IQ'^^) 

A— ►oo 

for all E D and A; G N. Due to the exponential boundedness of the Hermite functions, the 
moments {ip\Q^ip), A; G N, of the probability measure determine it uniquely [9]. Since D is 

a dense subspace, the probability measures p^, -0 G -D, determine, by polarization, the spectral 
measure Q of Q. To conclude that on the basis of the statistical data ([3]), the observable E^ 
would converge to Q, one needs to know that also E'^ is determined by its moment operators 
E^[k], k & N, on D. Again, for all G D, the probability measures are exponentially 



^Some of the technical details behinds these computations have been studied in |17j . 
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bounded, so that each is determined by its moments p^ [k] = {ilj\E^[k]ilj) , k G N. Hence, 
by polarization, is determined by the numbers {ip\E'^[k]ilj) , k & N,ip & D. 

Let now A„, n G N, be an increasing sequence of the coupling constants, with A„ — > oo, and 
let (E") „gN be the sequence of the semispectral measures E'^". The above results show that Q 
is the moment limit of the sequence {E"-)nm on D, that is, we may write 

(4) lim = Q 

n— >oo 

(on D in the sense of moment operators), for further technical details, see p4]). We remark 
that in this case also the effects E'^{X) tend weakly to the projections Q(X) for all X G B{M) 
whose boundaries X (1 X' are of Lebesgue measure zero, [T4] . 

The corresponding limits for the case Var {Qq, 0) ^ can similarly been worked out, for 

instance, if is chosen to be the Gaussian state 0„(x) = (^)^^^ e~"^^/^, and one considers the 
limit 77, — »• cx). 

3.1.2. Indirectly measured observable. In addition to obtaining the limit (J!]), formula ([2]) can 
also be solved directly for the numbers {ip\Q^ip) , G Z^, G N. Indeed, one may write recursively 

(5) (^igV) = {m[m - E k G N. 

These numbers are the moments of the probability distributions for p = Due 
to the exponential boundedness of these distributions they are uniquely determined by their 
moments (iplQ^ip), k eN, and by the density of D, the polarization identity then implies that 
this statistics is sufficient to determine also the position observable Q. Though the actually 
measured observable in this model is the unsharp position E'''''^, the measurement statistics 
allows one to determine also directly, without any limit considerations, the 'unobserved' sharp 
position Q. In Sect. [6] we discuss still another method to obtain the position distribution p^, 
p = \ip){ip\, from the actually measured distribution f * p^ hy inverting the convolution. 

3.2. The balanced homodyne detection observable. The balanced homodyne detection 
scheme is a basic measurement scheme in many quantum optical applications, including contin- 
uous variable quantum tomography as well as continuous variable quantum teleportation. Such 
a measurement scheme determines an observable E^ which depends on the coherent state \z), 
2; G C, of the auxiliary field. An important property of these observables is that on the level of 
statistical expectation values they agree with the quadrature observables Qe = ;^(e~*^a+e*^a*), 

z = re*^, of the relevant field mode, with the annihilation operator a. The explicit structure of 
these observables E^ has been studied in great detail and, in particular, their moment operators 
are determined [T4] . 

To express the relevant results here, we let -D(a) stand for the domain of the annihila- 
tion operator (which, in terms of the fixed number basis {|n)}„gN C H, is D{a) = {ip E 
^gpjTz |(n|95)P < 00}), and — a*a is the corresponding (selfadjoint) number operator. 
The first and the second moment operators of such a balanced homodyne detection observable 
E^ , z = re*^, are known to be as follows: 

-S^[l]|D(a) = Qe\D{a), 

E^[2]\n(a.) = (g,|^(,))2 + V^iV. 
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Here e.g. [A;] denotes the restriction of the /c*^ moment operator E^[k] of to the 
domain D[a'^), k = 1,2. By definition, the noise operator N{E^) has the domain D{N{E^)) = 
D{E[2]) n D{E[lf), which includes the set D{N) = D{a?) because of the above operator 
relations. Hence, \r~'^N C N{E^). But N{E^) is symmetric and selfadjoint, so that 
N{E^) = ^r~^N. This would again suggests that in the limit r = \z\ ^ oo, the intrinsic 
noise N[E^) goes to zero and thus the measured observable would approach the quadrature 
observable Qg. Like in the previous case. Sect. 13. ll this limit requires further considerations. 

Actually, the restrictions of all the moment operator E^[k] on the domains D{a^), /c G N, 
can be determined, and they are of the form 

(6) E'MDia^) = mDia^)f + ^C'fc(r,^^), 

where Cfc(r, 6) = J2n m <j. Cn,m(^5 6')(a*)"a'", and each is a bounded complex function on 
[1, oo) X [0, 27r) [14]. Let Dcoh = lin{|w) | w e C}, so that Dcoh is a dense subspace contained in 
all D{a^), A; G N. For each unit vector -0 G Dcoh, the probability measure p^"^ is exponentially 
bounded so that it is determined by its moment sequence ((il)\E^[k]'il)) ,k G N. Since Dcoh is 
dense, these probability measures define again the whole operator measure E^ [T4] . 

Let now (r„) be a sequence of positive numbers converging to infinity. For this choice, let 
Zn{(^) = TnC*^, where the phase 6 G [0, 27r) is also fixed, and let E^ be the corresponding balanced 
homodyne detection observable. By the above results it now follows that the spectral measure 
Qg is the only moment limit of the sequence of observables Moreover, for any unit vector 

ip, limn^oo P^" {X) = p^^lX) for all X G B{R) whose boundary X fl X' is of Lebesgue measure 
zero [14]. In this sense one can say that the high amplitude limit of the balanced homodyne 
detection scheme serves as an experimental implementation of a quadrature observable. 

Again, one may solve the statistical moments {ip \ Qgip) from ([6]) for all ip G D{a^). However, 
in this case they are not directly expressible in terms of actually measured moments {ip \ E^[k]ilj). 
The high amplitude limit is needed for that end. 

To close this section we mention that in a recent paper Man'ko et al |24j has proposed 
to use the first and second moments of the measurement statistics of the (limiting) balanced 
homodyne detection observables associated with the phases 6,6 + ^ + f , to empirically test 
the uncertainty relations for the conjugate quadratures (associated with 6,6 + ^). Clearly, for 
any ip G D{a^), with the choice 6 = and notations Qo = Q, Qil = P, 

Var^(EOVar^(E*0 = ((^12]) - {E'[l]f) ((^-[2]) - (i^^^l])') 

= {{Q') + \ r~HN) - (Q)2) ((P2) + i r-^N) - (P)') 

= (Var^Q) + ir-2(iV)) (Var^(P) +'^r~^{N)) > i, 

which allows one to test the statistics in this respect for any \z\ = r. The marginal statistics 
of the limiting eight-port homodyne detection observables of Section 14.31 leads to a similar 
inequality, except with the lower bound 1. We wish to point out that the test proposed in [24] 
is actually an experimental check for the correctness of the quantum mechanical description of 
balanced homodyne detection, since any violation of the above inequality would suggest that 
the description is incorrect. 
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4. Combining position and momentum measurements 



We shall go on to combine the above measurement schemes to produce sequential and joint 
measurements for position and momentum. We consider first the sequential application of a 
standard position measurement with any momentum measurement. Sections 14.21 and [43l deal 
with the Arthurs-Kelly model and the eight-port homodyne detection scheme. 

4.1. Sequential combination. Consider an approximate position measurement, described by 
the von Neumann model, followed by a sharp momentum measurement. This defines a unique 
sequential joint observable, a covariant phase space observable G^''^ : i3(R^) L(T-C), with the 
marginals 

= ixx*e)iQ), 
Gl^\Y) = {xY*f){P). 

Here we have the probability densities e{q) = A|0(— Ag)p and f{p) = — f )P, where (j) 
||</)|| = 1, is the initial probe state, and denotes the Fourier transform of 0. If G Cj^(R), we 
have G D{Qq) fl D{Pq) for each G N, in which case the moment operators of the marginal 
observables are 

(7) G'Ak] = X^(^)A-(0|Q^o</')Q'-^ 

j=0 ^ ^ 

(8) G'Ak] = f^V^(0|p<;0)p'=-^ 

As shown before, we have 

lim {^lj\G^,''^[k]^lj) = (V^IQV) 

A^oo 

for all G Cj^(]R). In the case of the second marginal we see that for any ip G Cj^(M) there are 

values of A; G N for which {ip\G2''^[k]ip) tends to infinity as A increases. For example, the limit 
of the second moment is never finite since (0|Po^0) is always non-zero. That is, the limits of the 
moments of the probability measure X h-^ {tlj\G2''^{X)tp) = {tlj\G^''^(R x X)-0) are not moments 
of any determinate probability measure, and hence they do not determine any observable. 

Another way to look at the limits of the marginal observables is to choose a sequence of 
initial probe states (0n)neN C ^^^(IR), such that |0„p approaches the delta distribution as n 
increases. For example, choose the Gaussian states 

0n(a;) =y-) e "2 , 

in which case the explicit forms of the moment operators and '"[^1 can easily be 

computed: 

(9) GJ-W ^ t f!);^r(i±i)Q- 



i=0, % even 



(10) GY\k\ 




where T denotes the gamma function. Taking the limit n ^ oo one gets a result similar to the 
one considered before (A oo). 

As expected, the limit procedures cannot give both the Q and P -distributions, but as it is 
obvious from ([THH]) and flOlfTOll the method of moments can again be used. We return to that in 
Sect. El 

Again, the convolution structure allows one to easily compute the distances between the 
marginals and the sharp position and momentum observables. One finds that 

d{Gl''^,Q) = jj \xmx)\^dx, 

d{G^''^,P) = ^ j \x\\^{x)\^dx, 

showing, that the product of the distances does not depend on A. Since the distances are Fourier- 
related, their product has a positive lower bound, that is, inf^^g^^^ d{G\''^ ^ Q) ■ d{G2''^, P) > 0. 
For example, in the case of the Gaussian initial states 0„ one has d{G{^ , Q) ■ (/(Gj P) = ^ fo'^ 
all n e N. 



4.2. Arthurs-Kelly model. The Arthurs-Kelly model [2J as developed further by Busch [4] 
(see also [28i, i29j) is based on the von Neumann model of an approximate measurement. It 
consists of standard position and momentum measurements performed simultaneously on the 
object system. Consider a measuring apparatus consisting of two probe systems, with associated 
Hilbert spaces Tii and 7^2- Let 0i ® 02 £ '^^i ® 'H2 be the initial state of the apparatus. The 
apparatus is coupled to the object system, originally in the state G 7i, by means of the 
coupling 

-iAQ(giPi (gi/2 pi/^PfSi-fl ®Q2 

which changes the initial state of the object-apparatus system \l/o = ^®0i ®02 into = U'^q. 
The final state has the position representation 

y, z) = ^{x + iJ^z)(j)i{y - \x)(p2{z). 

Notice, that the coupling ( ITTi ) is a slightly simplified version of the one used by Arthurs and 
Kelly. However, it does not change any of our conclusions. 

The measured covariant phase space observable G is determined from the condition 

{^\G{X X Y)^) = (^|/ ® Q(AX) ® PifiY)^), 

for all X,Y E i3(R), and the marginal observables Gi and G2 turn out to be 

(12) Gi(X) = (xx*(eo*l4"¥))(Q), 

(13) G2{Y) = (xy*(/o*|0S'^r))(^), 

where Cq and /o are the probability distributions related to the original single measurements, 
i.e. eo{q) = A|0i(— Ag)p and fo{p) = /i|02(— /^p)P, and we have used the scaled functions 
= "\/A0i(Ag) and (f)'f\p) = ■^02(^) • If we choose the initial state of the apparatus to 
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(11) 




Figure 1. The eight-port homodyne detector 



be such that 0i,</)2 € Cj^(R), the moment operators can be computed: 

(14) G,[k] = X^Xl(^)(^)A-("-)(-/x)^(0iiQr>i)(02ig^02)g 

^— n i—n V / V / 



n=0 i=0 

(15) G2[k] = X;XlHM/i-("-)(-A)^(02|Pr>2)(0i|P^0i)P'^^^ 



n=0 i=0 



It is clear from equations (fT2] -[T3l). that the Q- and P-distributions cannot be simultaneously 
obtained as limits of the marginals, since the distributions Cq * {(f'i^^l'^ ^ind /o * cannot 
both be arbitrarily sharply concentrated. However, equations (1 1411 151) show that the method of 
moments can be used. 



4.3. Eight-port homodyne detector. The eight-port homodyne detector [201 El] consists of 
the setup shown in Figure [H The detector involves four modes and the associated Hilbert spaces 
will be denoted by Hi, 7^2, '^s and H^. Mode 1 corresponds to the signal field, the input state 
for mode 2 serves as a parameter which determines the observable to be measured, and mode 
4 is the reference beam in a coherent state. The input for mode 3 is left empty, corresponding 
to the vacuum state. We fix a photon number basis G N} for each TCj, so that the 

annihilation operators aj, as well as the quadratures Qj = :^('2j + ctjOi Pj = ':)2^^*j ~ 
the photon number operators A'^ = a*aj are defined for each mode j = 1, 2, 3, 4. 

The photon detectors Dj are considered to be ideal, so that each detector measures the sharp 
photon number Nj. The phase shifter is represented by the unitary operator e*^^", where ^ is 
the shift. There are four 50-50-beam splitters B12, -B43, f/13, f/24, each of which is defined by 
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its acting in the coordinate representation: 



(16) L2(r2) ^^^J^ (a;, + x,),^(-Xi + x,))) e L\R^). 

In the picture, the dashed line in each beam splitter indicates the input port of the "primary 
mode", i.e. the mode associated with the first component of the tensor product L^(]R)(g)L^(]R) ~ 
L^(M^) in the description of equation (flGll . The beam splitters are indexed so that the first 
index indicates the primary mode. 

Let \V2z) be the coherent input state for mode 4. We detect the scaled number differences 
^\^i3 and j^^24) where = Ii®Nj — Ni^Ij, so that the joint detection statistics are 
described by the unique spectral measure extending the set function 

(X,r) ^ p\'\-'^i3{X) ® p\'\''^24(Y) = Di(X) ® D2(F), 

where the operator acts on the entire four-mode field. 

Let p = |^)(^/'| and a be the input states for mode 1 and 2, respectively. Then the state of 
the four-mode field after the combination of the beam splitters and the phase shiter is 

W^P,<x,.,e = Un ® f/24 {Buip ® ^)5*2 (g) \z){z\ (g) \ze'^){ze'^\ ) U*.^ ® U;^. 

We regard a, \\/2z) and ^ as fixed parameters, while p is the initial state of the object system, 
i.e. the signal field. The detection statistics then define an observable G^''^'^ : i3(M^) L{Hi) 
on the signal field via 

Tr[pG'^''^'«(X X Y)] = Tr[iy,,.,,,5Di(X) ® D,{Y)]. 

This is the signal observable measured by the detector. 

Let denote the covariant phase space observable generated by a positive trace one operator 
T, that is, 

(17) G^iZ) = ^ [ W,,TW;/qdp 

for all Z G yB(M^), where Wqp, {q,p) G M^, are the Weyl operators associated with the posi- 
tion and momentum operators Q and P. Let C : H2 ^ Ti-i denote the conjugation map, i.e. 
{Cip){x) = ip{x) in the coordinate representation, and let (r„) be any sequence of positive num- 
bers tending to infinity. It was shown in [15] that the measured observable G^""'^'2 approaches 
with increasing n the phase space observable generated by CaC~^, that is, 

lim G''"'"'5(Z) = G^''^"\Z) 

n— >oo 

in the weak operator topology, for any Z G B{R^) such that the boudary Z H Z' has zero 
Lebesque measure. 

In general, it is difficult to determine the domains of the moment operators of the covari- 
ant phase space observable G'-^"'-^ . However, if the generating operator CaC~^ is such that 
QWCaC-^ and PWCaC-^ are Hilbert-Schmidt operators for all k E N, then according to 

10 



[T6l Theorem 4] we have 
(18) 

(19) G^^'^<'-\k] 

5. Simultaneous measurements of Q and P 

In the three different measurement models considered above, the actually measured observ- 
able is a covariant phase space observable for an appropriate generating operator T. Hence, 
the marginal observables G\ and G'^ are convolutions of the sharp position and momentum 
observables with the Fourier related probability densities / and g defined by T, respectively. 
Indeed, if T = ^j^j|^7j)(?7i| is the spectral decomposition of T, then /(g) = X]i^«l'7i(~5')P 
dip) = J2i'^i\Vi{~P)\'^ ■ Due to this structure, the moment operators of the marginal observables 
and G^ can be written in simple forms as polynomials of either Q or P. That is, for any 

k 

i=0 
k 

1=0 

where the coefficents and sf^ depend on the model in question and s^q = s^q = 1 in each 
case. From these, the recursion formulae for the moments of the position and momentum 
distributions p'^ and p^, with p = of the object to be measured can be computed: 

k 

(20) (^igV) = ii^iGim-Y^sZm'-'i^), 

1=1 

k 

(21) (v^IpV) = {^p\Glm-J2sM{i'\P'-'i^)- 

1=1 

If ip is chosen to be, for example, a linear combination of Hermite functions, the distributions 
and p^ are exponentially bounded and as such, are uniquely determined by their respective 
moment sequences {{ip\Q^4'))km and {{ip\P^ip))k(^N- In this sense one is able to measure simul- 
taneously the position and momentum observables Q and P in such a vector state in any of the 
three single measurement schemes collecting the relevant marginal information. Furthermore, 
since the linear combinations of Hermite functions are dense in L^(M), their associated distri- 
butions and p^ suffice to determine the whole position and momentum observables Q and 
P as spectral measures. 

6. Concluding remarks 

We have shown with three different measurement models that the statistical method of 
moments allows one to determine with a single measurement scheme both the position and 
momentum distributions p^ and p^ from the actually measured statistics for a large class of 
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E(^)(-irTr[aQ^]Qt-" 

n=0 ^ ^ 
n— n \ / 



initial states p. In each case the actually measured observable is a covariant phase space 
observable G'^ whose generating operator T depends on the used measurement scheme. Such 
an observable is known to be informationally complete if the operator T satisfies the condition 
Tr[PVgpT] 7^ for almost all (g,^) G ]R^[lJ. Recently it has been shown that this condition 
is also necessary for the informational completeness of [THj. Neither the used models nor 
the method of moments depend on this assumption. Indeed, if, for instance T = \ri){ri\, with a 
compactly supported 77, so that G^ is informationally incomplete, the equations (l20l-[2Til can still 
be used to determine and provided that these distributions are exponentially bounded, for 
instance if p = with ip in the linear hull of the Hermite functions. If, however, the phase 

space observable G^ is informationally complete and if one is able to reconstruct the state p 
from this informationally complete statistics 1i[pG'^ {Z)], Z e S(M^), then, of course, one knows 
the distribution of any observable, in particular, the position and momentum distributions 
and p^ . However, the reconstruction of the state from such a statistics is typically a highly 
difficult task, see e.g. [27]. In the special case of the generating operator T being the Gaussian 
(vacuum) state T = |0)(0|, the distribution Z ^ Tr[pGl°^(Z)] is the Husimi distribution of 
the state p. For that, a reconstruction formula is well known and simple [^. Indeed, writing 
z = j^iq + w). one has W,p\0) = \z), and Tr[pG'|o)(Z)] = J^Qpiz)d^z, with Qp{z) = l{z\p\z) 
being the Husimi Q-function of the state p. Using the polar coordinates, the matrix elements 
of p with respect to the number basis are 

~ {2n + k)\ rfr2"+fc ' 

where 

'•277 



It is to be emphasized that the reconstruction of the state requires, however, full statistics of 
the observable C'''^ The marginal information, which is used in the method of moments, is 
clearly not enough to reconstruct the state even in the case where the position and momentum 
distributions are exponentially bounded. To illustrate this fact, let us consider the functions 
^a,b{q) = (^)^'''^e"("+*^)''^ with a,6 G M, a > 0. The Fourier transform of ipa,b is 

and the position and momentum distributions are 

\fa,b{q)\ 



2a ^ 



\n J 

which are clearly exponentially bounded. For 6 7^ 0, we see that pi = \fa,b){'^a,b\ and p2 
\ipa-b){'^a-b\ are different states, but pf = p^ and pf = pf. The marginal probabilities are 

pfiY) = lyig * pf)iy)dy = Iy{9 * P^){y)dy =pf{Y), 

12 



for all X,Y E B{R), with g{x) = -^e ^ , so the marginal distributions are equal. It follows 
that the state cannot be uniquely determined from the marginal information only. 

Since the marginal observables Gj and are of the convolution form with densities, the 
position and momentum distributions can also be obtained if one is able to invert the con- 

qT qT 

volution. Indeed, for any initial state p = \ip){ip\ the marginal distributions pp^ and 
have the densities f * p'^ and g * p^ , where /(g) = g)P, g{p) = Sj^il^j(— P)P) with 

T = J2i'^i\Vi){Vi\j ^iid p*^ = \ip\'^, p^ = l^/'p. The unknown distributions p*^ and p^ can be 
solved from the measured distributions f *p^ and g* p^ by using either the Fourier inversion or 
the differential inversion method. Like the method of moments, these methods have their own 

specific restrictions. In fact, by the Fourier theory, one has, for instance, / * = v^27r/-p'3, so 

that p'3 = (27r)^^/^/ * p*^//, provided that / is pointwise nonzero. Iff* p^ / f is an L^-function, 
then the function 

1 . — - — , 

— / e^'f*pQ{t)ij{t)dt 

coincides with the distribution p^ (almost everywhere). Obviously, this puts strong restrictions 
on the actually measured distribution j * p^ as well as on the 'detector' density / = /(T). 
The method of differential inversion is known to be applicable whenever the detector densities 
/ and g have finite moments [12]. In the special case of T = |0)(0|, so that / and g are the 
Gaussian -^e~^ , one has 

fc=0 

k=0 ^ 

provided that the right hand sides exist [12], which is a further condition on the initial state p. 

To conclude, the statistical method of moments provides an operationally feasible method to 
measure with a single measurement scheme both the position and momentum distributions p^ 
and p^ for a large class of initial states p, the relevant condition being the exponential bound- 
edness of the involved distributions. This method requires neither the state reconstruction nor 
inverting convolutions. 



Appendix A. Proof of lemma [T] 

If P is a projection in the range of M, then P commutes with any effect M{Z), Z e i3(M^) 
(see, for instance, [22l Th. 1.3.1, p. 91]). Therefore, the marginals Mi and M2 are mutually 
commutative, i.e. Mi(X)M2(F) = M2{Y)Mi{X) for all X,Y e B{R), and the map {X,Y) ^ 
Mi{X)M2{Y) is a positive operator bimeasure, and extends uniquely to a semispectral measure 
G : i3(R2) L{n), with G{X x Y) = Mi{X)M2{Y) for all X,Y e B{R) (see, e.g. , Theorem 
1.10, p. 24, of [3j). Let X,Y e i3(M). Since Mi(X) and M2{Y) commute and one of them is a 
projection, we have G{X xY) = Mi{X)M2{Y) = Mi{X) A M2{Y), the greates lower bound of 
Mi(X) and M2{Y), [251 Corollary 2.3]. Since also M{X x F) is a lower bound for Mi(X) and 
M2IY), we obtain M{X xY) < G{X xY). It follows that M{Z) < G{Z) for any Z eJ^, where 
J-' is the algebra of all finite unions of mutually disjoint sets of the form X x Y , X,Y E B{R). 
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Denote M = {Z e B{R^) \ M{Z) < G{Z)}. Now M is a monotone class. [If is an 
increasing sequence of sets of A4, then for any 9? G Ti, we have 

{ip\M{UnBM - {ip\G{UnBM = ]im{{ip\M{BM ' {^\G{BM) < 

n 

because e.g. Z 1— > {(f\M{Z)(fi) is a positive measure. This shows that U„S„ e A4. Similarly, we 
verify the corresponding statement involving decreasing sequences, and thereby conclude that 
is a monotone class.] Since T G /A, and T is an algebra which generates the cr-algebra 
BiR"^), it follows from the monotone class theorem that M{Z) < G{Z) for all Z e B{R'^). Let 
Z e B(M^), and let (p & Tl be any unit vector. Since M^^^ and G^^ip are probability measures, 
we get 

1 - M^,^(Z) = M^,^(R2 \ z) < G^,^{R^ \Z)^1- G^,^{Z), 

implying that {(p\G{Z)ip) < {(fi\M{Z)(f) . Since (p was arbitrary, this implies G{Z) < M{Z). 
The proof is complete. 

□ 



References 

[1] S.T. Ali, E. Progovecki, Classical and quantum statistical mechanics in a common Liouville space, Physics 
89A (1977) 501-521. 

[2] E. Arthurs, J. Kelly, On the simultaneous measurements of a pair of conjugate observables. Bell System 
Tech. J. 44 (1965) 725. 

[3] C. Berg, J.P.R. Cliristonscn, P. Rcssel, Harmonic Analysis on Semigroups, Springer, Berlin, 1984. 

[4] P. Busch, Unbestimmtheitsrelation und simultane Messungen in der Quantentheorie, Ph.D. thesis. Uni- 
versity of Cologne, 1982. Enghsh translation: Indeterminacy relations and simultaneous measurements in 
quantum theory, Int. .J. Theor. Phys. 24 63-92 (1985). 

[5] P. Busch, M. Grabowski, P. J. Lahti, Operational Quantum Physics , Springer, Berlin, 1995. 

[6] P. Busch, J. Kukas, P. Lahti, Measuring position and momentum together. Physics Letters A 372 (2008) 
4379-4380. 

[7] G. M. D'Ariano, C. Macchiavello, M. G. A. Paris, Detection of the density matrix through optical homodyne 

tomography without filtered back projection, Phys. Rev. A 50 (1994) 4298. 
[8] D.A. Dubin, J. Kiukas, .J. -P. Pellonpaa, private communication 2008. 
[9] G. Freud, Orthogonal Polynomials, Akademia Kiado, Budabest, 1971. 
[10] W. Heisenberg, Uber den anschauHchen inhalt der quantentheorischen kinematik und mechanik, Z. Phys. 
43 (1927) 172-198. 

[11] I. I. Hirschman, D. V. Widder, The Convolution Transform, Princeton University Press, Princeton, 1955. 
[12] R. G. Hohlfeld, J. I. F. King, T. W. Drueding, G. V. Sandri, Solution of convolution integral equations by 

the method of differential inversion, SIAM J. Appl. Math., 53 (1993) 154-167. 
[13] A. S. Holevo, Covariant measurements and uncertainty relations, Rep. Math. Phys. 16 (1979) 385-400. 
[14] J. Kiukas, P. Lahti, On the moment hmit of quantum observables, with an application to the balanced 

homodyne detection, J. Mod. Optics 55 (2008) 1175-1198. 
[15] J. Kiukas, P. Lahti, A note on the measurement of phase space observables with an eight-port homodyne 

detector, J. Mod. Optics (2007) 
[16] J. Kiukas, P. Lahti, K. Ylinen, Phase space quantization and the operator moment problem, J. Math. Phys. 

47 (2006) 072104/18. 

[17] J. Kiukas, P. Lahti, K. Yhnen, Semispectral measures as convolutions and their moment operators, J. 

Math. Phys. 49 (2008) 112103/6. 
[18] J. Kiukas, R. Werner, private communication 2008. 

14 



[19] P. Lahti, J. -P. Pellonpaa, K. Ylinen, Two questions on quantum probability, Phys. Lett. A 339 (2005) 
18-22. 

[20] U. Leonhardt, H. Paul, Phase measurement and Q function, Phys. Rev. A 47 (1993) R2460-R2463. 
[21] U. Leonhardt, Measuring the Quantum State of Light, Cambridge University Press, Cambridge, 1997. 
[22] G. Ludwig, Foundations of Quantum Mechanics 7, Springer- Verlag, Berlin 1983. 

[23] A. Luczack, Instruments on von Neumann algebras, Institute of Mathematics, Lodz University, Poland, 
1986. 

[24] V.L Man'ko, G. Marmo, A. Simoni, F. Ventriglia, A possible experimental check of the uncertainty relations 

by means of homodyne measuring photon quadrature . iiarXi v :0811.4115i yl. 
[25] T. Moreland, S. Gudder, Infima of Hilbert space effects. Linear Algebra and its Applications 286 1-17 

(1999). 

[26] M. Ozawa, Quantum measuring processes of continous observables, J. Math. Phys. 25 (1984) 79-87. 
[27] M. Paris, J. Rehacek (Eds), Quantum State Estimation, Lect. Notes Phys. 649, Springer- Verlag, Berlin, 
2004. 

[28] M. G. Raymer, Uncertainty principle for joint measurement of noncommuting variables. Am. J. Phys. 62 

(1994) 986-993. 

[29] P. Torma, S. Stenholm, L Jex, Measurement and preparation using two probe modes, Phys. Rev. A 52 

(1995) 4812-4822. 

[30] J. von Neumann, Mathematische Grundlagen der Quantenmechanik, Springer- Verlag, Berlin,, 1932. 
[31] R. Werner, Quantum harmonic analysis on phase space, J. Math. Phys. 25 (1984) 1404-1411 . 
[32] R. Werner, Dilations of symmetric operators shifted by a unitary group, J. Func. Anal. 92 (1990) 166-176. 
[33] R. Werner, The uncertainty relation for joint measurement of position and momentum, Qu. Inf. Comp. 4 
(2004) 546-562. 

Institute for Mathematical Physics, Technical University of Braunschweig, Braunschweig, 
Germany 

E-mail address: jukka.kiukas@utu.fi 

Department of Physics and Astronomy, University of Turku, Turku, Finland 
E-mail address: pekka.lahti@utu.fi 

Department of Physics and Astronomy, University of Turku, Turku, Finland 
E-mail address: jussi.schultzOutu.f i 



15 



